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WEATHER PREDICTION METHOD 
FOR FORECASTING SELECTED EVENTS 

STATEMENT REGARDING FEDERALLY SPONSORED 
5 RESEARCH OR DEVELOPMENT 

[0001] This invention was made in the performance of work by an inventor 
employed by the U. S. Government. The United States Government may have 
certain rights to this invention and a right to use any software products 
10 developed thereunder. Under authority of 35 U.S.C. §205, Federal agencies 
may withhold any information pertaining to this invention from the public for a 
reasonable time in order for a patent application to be filed. Furthermore, 
Federal agencies shall not be required to release copies of any software that is 
part of this invention. 

15 

CROSS-REFERENCES TO RELATED APPLICATIONS 

[0002] The present patent application claims the benefit of the filing date of 
U.S. Provisional Patent Application No. 60/413,936, filed September 26, 2002, 
20 and entitled "Weather Prediction Method for Forecasting Selected Events". 

FIELD OF THE INVENTION 

[0003] The present invention relates to the area of weather forecasting and 
25 modeling, and more particularly, to methods, systems, and computer program 
products for a weather prediction system that will produce short term probability 
forecasts for the occurrence of weather events of interest. 
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BACKGROUND OF THE INVENTION 

[0004] The weather has a profound effect on the personal and economic 
lives of people. Its pervasiveness affects various decisions made by individuals 
5 on how to order their affairs, optimize economic activity, and plan for future 
events. Different weather events have different significance, depending upon 
the activity involved. For example, lightning is a weather event that is known to 
be responsible for deaths and injuries, forest fires, personal and business 
property losses, and airline delays. The ability to predict the probability of 

10 lightning strikes is of great importance since it allows informed decisions to be 
made and proactive measures to be taken, in order to mitigate the negative 
economic and human impact resulting from these events. For example, forest 
fire fighters might be dispatched to areas with high probability of lightning strikes 
so that they could rapidly handle any fire that might be caused by the lightning 

15 strikes. The same scenario would hold true for the prepositioning of repair 
crews and equipment for electrical power outages caused by lightning strikes. 
Airlines would be able to better schedule and reroute flights to accommodate 
adverse circumstances caused by lightning strikes. 

[0005] Similarly, all other weather-related events, such as snowfall, heavy 
20 rain and flooding, large hail, damaging wind, and tornadoes, may have different 
but equally significant affects on various activities, human and economic. Many 
resources have been devoted to improving the ways in which the probability of 
occurrence of these weather-related events can be forecast. In order to 
forecast such weather-related events, data must be collected and analyzed on a 
25 real-time basis. The data collection may produce an enormous database that 
may exhibit little uniformity as to physical spacing, accuracy, and time 
sequencing of the measurements. The accuracy and extent of these data may 
vary from region to region depending upon the facilities available and the 
importance of having an accurate forecast to those personnel in the region. 
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[0006] Data must be obtained from which to forecast the probability of 
weather events for a given region. Currently, the Rapid Update Cycle 2 (RUC2) 
model is one method used for measuring the atmosphere across the United 
States and providing a useful data set for prediction of weather events. 
5 Becoming operational in 1998, the RUC2 model represented the first time that 
an hourly, 3-dimensional data set became available. It contains the parameters 
necessary for evaluating the potential for convection as well as heavy snow, 
rain, severe weather or any other selected weather event, i.e. pressure, 
temperature, moisture, and wind parameters. Another such model is the ETA 

10 Model, one of the operational numerical forecast models run at the National 
Centers for Environmental Prediction (NCEP). The ETA model is run four times 
daily, with forecast output out to 84 hours. Yet another model under 
development is the Weather Research and Forecasting (WRF) model. 
[0007] There have been numerous methods developed for predicting 

1 5 weather-related events. One recent effort at lightning prediction for the zero to 
three-hour period is the System for Convection Analysis and Nowcasting 
(SCAN; Smith et al. 1998). This technique is an extrapolative-statistical method 
based on three years of lightning, radar, and satellite developmental data. The 
algorithms use extrapolated forecasts of radar reflectivity, satellite-estimated 

20 cloud-top temperatures, and current lightning strike rates. However, these 
algorithms do not extend very far into the future because they rely so heavily on 
the extrapolation of current conditions. They also are unable to predict when 
and where new storms will develop. With the exception of these extrapolative 
forecasts, no short-term (0 to 6 hour) lightning forecasts are in operational use. 

25 In addition, previous attempts, using Model Output Statistics or MOS, have tied 
predictive equations to specific models, which cannot change as the models 
themselves change and improve. Other similar methods are given according to 
U.S. Pat. No.6,405,134, to Smith, entitled "Method and Apparatus for Predicting 
Lightning Threats Based on Radar and Temperature Data", that provides 

30 forecasts of one hour or less. 
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[0008] While sophisticated computer models predict many surface 
parameters, such as surface wind (10 meter), surface temperature, and dew 
point (2 meter), other important quantities are not forecast directly. Fields 
such as amount of cloud cover; ceiling; visibility; probability and type of 
5 precipitation; thunderstorm probability; and severe thunderstorm probability 
are not explicitly predicted by most models. Two basic methodologies have 
been used to compute fields not specifically forecast by a model. These are 
known in the art as the Model Output Statistics (MOS) and the Perfect Prog 
(prognosis) methodologies. 

10 [0009] According to the Perfect Prog methodology, a statistical 
relationship is developed between the variable to be estimated and selected 
variables which can be forecast by a model. The predictors and predictand 
are observed quantities in the developmental sample. This relationship is 
applied using the output fields from the forecast models as predictors to 

15 estimate the predictand, treating the model forecast as a perfect prognosis. 
As improvements to the model are made, forecasts of the predictors should 
improve with corresponding improvement of the Perfect Prog forecasts. 
[0010] According to the MOS methodology, a statistical relationship is 
determined between the predictand and the predictors from a numerical 

20 model at a predetermined projection time. MOS equations in the past have 
been able to out-perform a Perfect Prog approach because they take into 
account model biases and the decline in the model accuracy with increasing 
forecast projection. The most serious drawback to the MOS methodology, 
however, is that it requires a stable numerical model. Given the rate with 

25 which changes and improvements are made to today's models, it is doubtful 
that any model will remain static long enough to develop new MOS 
equations. Alternatively, the use of the MOS methodology requires the 
newer version of any model to be "rerun" on a minimum of one year (and 
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usually two to three years) of model output data, which may require 
expenditure of extensive data collection and computational resources. 
[0011] Because today's complex computer forecast models change more 
frequently than those of the past, recent research in forecasting weather- 
5 related events has not been transferred to the MOS statistical forecast 
methodology. Other MOS-based statistical forecasts have ended when the 
particular model they were based upon was discontinued (e.g., the Limited 
area Fine Mesh (LFM) model). MOS equations derived from the older 
models such as the Nested Grid Model (NGM) have not been updated since 

1 0 the NGM code was "frozen" in 1 991 . 

[0012] The basic problem is one of data compression. In other words, vast 
amounts of historical and real time data are now available through various 
models. So much data is available at different levels and times that it is often 
difficult to reduce and/or compress the data in a way that is meaningful for the 

15 particular weather-related event of interest. What is needed is an improved 
statistical prediction system for more accurately predicting the occurrence of 
selected weather events; for reducing or compressing the data by identifying 
meaningful parameters that are pertinent to the weather-related event being 
examined; for producing forecasts that bridge the existing gap between 

20 extrapolative systems and model-based systems by using both analyses and 
model forecasts; and for improving the general understanding of environmental 
characteristics which support the significant weather events 

SUMMARY OF THE INVENTION 

25 

[0013] An object of the present invention is to provide a prediction system 
that will rapidly produce probability forecasts of selected weather events of 
interest from a set of predictors that account for a significant fraction of the 
known variance. 
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[0014] Another object of the present invention is to provide a system that is 
able to use input data provided by either Perfect Prog or MOS methodologies. 
[0015] Yet another object of the present invention is to develop a reduced 
set of uncorrelated predictors from the available predictor set, for a given 
5 geographical area, that more accurately forecasts selected weather events than 
heretofore. 

[0016] Yet another object of the present invention is to develop predictive 
equations according to the Perfect Prog methodology that can be applied to the 
ETA model, the RUC2 model, the Weather Research and Forecasting (WRF) 
1 0 model, or other forecast models. 

[0017] Yet another object of the present invention is to develop a weather 
variable prediction system that will provide forecasts for any selected type of 
weather event. 

[0018] Yet another object of the present invention is to provide a prediction 
15 system that will produce three-hour grid point probability forecasts of cloud-to- 
ground lightning, where the probability forecast assigns a probability of whether 
or not the lightning flashes occurring per hour is greater than a selected 
threshold value. 

[0019] Yet another object of the present invention is to provide a set of 
20 lightning climatological predictors that contain sufficient detail to represent for a 
selected weather-related event, the diurnal trends as well as intra-seasonal 
trends. 

[0020] Yet another object of the present invention is, for a given 
geographical area, to develop a set of predictive equations for lightning that can 
25 be applied to either observed gridded data (i.e. model analyses) or forecast 
data. 

[0021] To achieve the foregoing objects, and in accordance with the purpose 
of the invention as broadly described herein, the present invention provides 
methods, systems, and computer program products for a weather prediction 
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system that will produce probability forecasts of at least an hour for selected 
weather events. This method preferably comprises the following steps (not 
necessarily in the given order): defining a predetermined geographical area in 
terms of a regular grid, developing a data set of predictors consisting of a time 
5 series of measurements of a pre-defined collection of physical attributes, 
correlating the predictors to the predictands (to find the strongest correlation to 
predictands among the variables), performing Principal Component Analysis 
(PCA) to determine a much smaller, distinctive physical set of predictors to 
serve as new predictors for regression, and performing a logistic regression to 

1 0 determine the predictive equations for lightning occurrence. 

[0022] In an alternative embodiment, the steps of defining the geographical 
area in terms of a linear grid and developing a data set consisting of a time 
series of measurements of a pre-defined collection of physical attributes may be 
combined into the single step of inputting a data set obtained from any 

1 5 calculated gridded analysis of the data set and/or a model forecast. 

[0023] In another alternative embodiment, the system may be used to 
forecast the occurrence of any selected type of weather event, for example, the 
probability of heavy snow (such as 18 inches or more), the probability of fog 
(such as 1/4 mile visibility or less), probability of hail (such as 2 inches in 

20 diameter or greater), or the probability of a tornado (such as strength F2 or 
greater). The range of possible weather events that can be forecast by the 
system is only limited by the types of physical measurements of weather-related 
data that can be made that would reasonably be related to the weather event in 
question. 

25 [0024] The present invention will now be described with reference to the 
following drawings, in which like reference numbers denote the same element 
throughout. These and other features, aspects, and advantages of the present 
invention will become better understood with reference to the following 
drawings, description, and claims. These objects should be construed to be 



Page 7 of 48 



PATENT 
03-00076-03 



merely illustrative of some of the more prominent features and applications of 
the present invention. Many other beneficial results can be attained by applying 
the disclosed invention in a different manner or by modifying the invention within 
the scope of the disclosure. 

5 

BRIEF DESCRIPTION OF THE DRAWINGS 

[0025] FIG. 1 is a diagram showing in detail the method of calculating a 
forecast probability for a selected weather-related event using the inventive 
1 0 method, according to an embodiment of the invention. 

[0026] FIG. 2 is a plot of eigenvalues computed from the set of rotated 
Principal Components ordered by descending value. 

[0027] FIG. 3 is a diagram illustrating the approximate vertical level 
associated with different groups of data, according to an illustrative embodiment 

1 5 of the present invention. 

[0028] FIG. 4 is a flowchart of showing how the general weather prediction 
method may be applied when the weather-related event is lightning strikes, 
according to preferred embodiments of the present invention. 
[0029] FIG. 5 is a diagram showing the ten RPC terms for lightning 

20 prediction and the percentage of the variance explained by each term, where 
the total variance from all predictors contained in the ten RPCs is 59%, 
according to preferred embodiments of the present invention. 
[0030] FIG. 6 is a diagram showing the response of a binary predictand f(W) 
calculated from a logistic regression W, according to an embodiment of the 

25 invention. 
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DETAILED DESCRIPTION OF THE INVENTION 

[0031] The following detailed description shows the best currently 
contemplated modes of carrying out the invention. The description is not to be 
5 taken in a limiting sense, but is made for the purpose of illustrating the general 
principles of the invention and the best mode for practicing the invention, since 
the scope of the invention is best defined by the appended claims. 
[0032] The present invention defines methods, systems, and computer 
program products for a weather prediction system that will produce near-term 

10 probability forecasts of selected weather events for a specified area of interest. 
The invention employs the Perfect Prog methodology since it can be applied 
to different computer models at any forecast (or analysis) time and can be 
applied even as the models change. This advantageous feature makes the 
invention generically applicable to a large number of models. 

15 [0033] The invention employs a multivariate statistical method known as 
Principal Components Analysis, or PCA, to essentially compress the data so 
that it may be more easily worked with. The atmospheric sciences have 
employed PCA in the prior art to reduce data sets containing large numbers of 
variables to a smaller set of variables called principal components (PC), which 

20 still represents a large fraction of the variability contained in the original data 
set. However, PCA can achieve this reduction only if there is substantial 
correlation, or covariance, between the original variables. Each of the PCs has 
the property as a group of being uncorrelated. 

[0034] In addition to data compression, the PCs can also provide a better 
25 physical interpretation of atmospheric processes and can yield substantial 
insights into both the spatial and temporal variations of the fields being 
analyzed. However, each PC may represent a combination of several 
processes. In order to relate the PCs back to measurable variables in the 
original data set, the invention provides an inventive procedure of applying a 
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linear transformation, or "rotation", to the reduced set of PCs to provide a 
corresponding set of rotated PCs, or RPCs. Even though the original variables 
in the observed data set may be correlated, these RPCs retain all of the 
variability in the reduced set of unrotated PCs, but are uncorrelated if an 
5 orthogonal rotation is used; uncorrelated variables are a normal result of such a 
coordinate transformation. Often in multiple regression, predictor variables 
exhibit very large mutual correlations. The result is that the forecasting equation 
may perform very badly with this method. Because the PCs are uncorrelated, 
there is no multi collinearity. Therefore, the regression coefficients will be 

10 estimated with much better precision. 

[0035] While the PCA is extremely powerful for data compression and can 
yield substantial insights into both the spatial and temporal variation of the fields 
being analyzed, the linear transformation ("rotation") may be applied to assure 
the PCs correspond to a single physical process and do not reflect mixed 

15 processes. One of the properties of the PCA method is that the eigenvectors of 
symmetric matrices are perpendicular (orthogonal) to each other. According to 
the PCA method, the eigenvectors calculated from a symmetric data matrix 
serve to define an alternative coordinate system to view the data. However, the 
orthogonality property of the eigenvectors can cause problems with their 

20 interpretation. This new (rotated) coordinate system may be oriented such that 
each new axis may be aligned along the direction of the maximum joint 
variability of the data, consistent with that axis being orthogonal to the preceding 
axes. The first eigenvector may be the vector associated with the largest 
eigenvalue and is aligned parallel to the direction in which the data jointly exhibit 

25 maximum variation. While the first PC may represent an important atmospheric 
process, the maximal variance property of the PCA method makes it likely that 
other distinct atmospheric processes are also combined and mixed in with the 
first principal component or distributed across a number of other PCs. 
[0036] Since one purpose of the PCA method is to gain a better physical 

30 interpretation of the atmospheric processes leading to the atmospheric 
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phenomenon being examined, e.g. thunderstorms in the Western United States, 
it is desirable to apply the linear transformation (i.e. a varimax orthogonal 
rotation) to a subset of the original PCs, in order to form a new set of rotated 
PCs (RPCs). An advantage of applying the varimax orthogonal rotation is that 
5 the RPCs are physically consistent with subsets of the original data variables 
that coincide with specific atmospheric processes. In other words, the RPCs 
support the observed correlation between the original variables. Hence, the 
RPCs may provide a better understanding of the physical basis for the 
atmospheric process being studied. 

1 0 [0037] Suppose, for example, we were to assume (perhaps erroneously) that 
we had devised "the perfect set of predictors", where each predictor 
represented a distinct physical process in the atmosphere and each of the 
predictors was uncorrelated and independent. The PCA method might 
independently confirm or disprove the original predictor assumptions. Even if 

15 the underlying physical processes leading to the phenomenon being studied 
(such as how thunderstorms form) was not known or understood in advance, 
the RPCA method would be able to take a large set of predictors and assign 
them to a much smaller set of rotated vectors that could best describe the 
elemental atmospheric variables. 

20 [0038] The statistical calculations performed herein in support of the 
invention may be provided by any suitable statistical software package having 
PCA and logistic regression capabilities. The PCA techniques used in 
developing the results for the invention may be found in statistical software 
packages such as S-Plus 2000, which is available from Insightful Corporation 

25 (previously, MathSoft, Inc.), 1700 Westlake Avenue North, Seattle, WA) (1999). 
Another statistical software package containing routines useful for calculating 
portions of the inventive method may be the GEMPAK software, which was 
developed under government contract for use by the community of statistical 
users. 
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[0039] Referring now to FIG. 1, a detailed procedure for the PCA method 
may now be described. Briefly, a set of parameters related to the weather- 
related event of interest is identified, according to the block labeled 110. A set 
of observations of these parameters is made at a set of geographical points 
5 over a period of time, according to the block labeled 120. An n x m data matrix 
may be developed from the n observations of the atmospheric variable of 
interest for a particular area, where the data matrix contains m data parameters 
and this data matrix is scaled, according to the block labeled 130. 
[0040] A correlation matrix may then be computed from the scaled data 
10 matrix, according to the block labeled 140. The diagonal elements of the 
correlation matrix are the autocorrelations of each of the data parameters with 
themselves and are always equal to 1 . Since correlation r & between the /' and j 

variables is the same as the correlation r J( , then the correlation matrix will 

always be symmetric. The correlation matrix may be used rather than the 
15 covariance matrix, because the correlation matrix uses variance-normalized 
parameters, which assures that parameters with large differences in magnitude, 
such as potential vorticity (lO^roV'tffcg" 1 ) and the 500 mb height field (10 3 
meters), receive equal weight. 

[0041] According to the block labeled 150, the eigvenvalues (scalar 
20 quantities) and eigenvectors (a set of nonzero vectors) may then be derived 
from the correlation matrix. Next, the individual elements, also called 
coefficients or PC loadings, of each of the eigenvectors may be computed and 
those eigenvectors contributing to the most variance may be retained, according 
to the block labeled 160. Mathematically, there are as many eigenvectors as 
25 there are data parameters. However, if all the eigenvalues ordered by 
descending magnitude, assigned an ascending index number, and plotted, then 
only the eigenvalues that retain the largest amount of variance from the original 
data need to be retained. These would be those values where the curve has 
the steepest slope. For example, as can be seen from the illistrative graph 200 
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in FIG. 2, the first 10 PCs 210 contain the most information. The first 10 PC 
loadings 220 are the most significant when viewed in light of the rate of change 
of the slope on FIG. 2, and they may constitute the reduced set of PCs of 
interest. This reduced set of PCs may then be rotated by applying a linear 
5 transformation to give rise to a new set of vectors called the rotated PCs 
(RPCs), according to the block labeled 170. Each of these RPCs contains all of 
the original quantities, and each vector can be defined by the largest absolute 
value coefficients. 

[0042] The detailed steps for applying a PCA to the predictands are given in 
10 FIG. 1. The PCA may be initiated with the computation of the scaled data 
matrix [z] (according to the block labeled 130) and the correlation matrix [r] 
(according to the block labeled 140) from the data matrix 
[x]={a;.. :/ = l,...,/i;y = l,...,m}, where n is the number of observations and m is 

the total number of parameters. The scaled data matrix [z] may be computed 
15 from the original data matrix [x]according to the following equation: 

[z]=[z>]-' (wAm 

n 

where -[ljx] is the mean of [x], [l] is the matrix whose elements are all equal 

n 

20 to unity, and [Z)]" 1 is the diagonal matrix whose elements are the reciprocals of 
the sample standard deviations of the m parameters on the diagonal of [d]. 
The correlation matrix [r] can be computed from the scaled data matrix [z] by 
the following equation: 

25 M--4OTM 

rt-1 
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10 



[0043] Correlation matrix [r] may be used rather than the variance- 
covariance matrix because of the large difference in the magnitudes of the m 
parameters. If the data were not scaled, that is, normalized to lie between 1 
and -1, then the terms with the smaller variance might be incorrectly excluded. 
For example, because the variance of a set of 500 mb heights would normally 
be many times larger than the variance of a set of vertical velocity values, then 
the PCA covariance matrix might be dominated by the 500 mb heights. 
[0044] The scaled data matrix [z] may be related to the principal component 
scores [f] and the principal component loadings [a] by the following equation: 

[zMfM 



Later, this equation may be rewritten to solve for the principal component 
scores [f]. These scores may serve as the predictors in the regression to be 

1 5 presently described, in place of the original m parameters. 

[0045] In order to calculate the principal component scores, the principal 
component loadings [f] must first be computed. The principal component 
loadings, or weights, [a] may be computed according to the block labeled 150 
from the eigenvector matrix [e] and the diagonal matrix [a], the elements of 

20 which are the descending eigenvalues, A of the correlation matrix [r], as 
follows: 

25 The eigenvectors [e] and eigenvalues [a] may be calculated from the 
symmetric correlation matrix [r]. 
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[0046] Next, the eigenvalues for all m parameters are examined according to 
the block labeled 160 to determine which parameters contribute the most 
variance. One method making this determination is to order all the eigenvalues 
by descending magnitude, plotting them on a graph, and using an empirical 
5 criterion for dividing the parameters into those that are significant and those that 
are not. One criterion might be to discard the remaining parameters when the 
difference between two successive eigenvalues is less than a given amount. 
Another criterion might be to examine the slope of the curve, as shown in the 
illustrative graph 200 of FIG. 2, and discard parameters to the right of the point 
10 where the slope approaches zero. At any rate, the PCs corresponding to the 
eigenvalues representing the greatest variance may be chosen for the reduced 
set of PCs. 

[0047] The reduced set of p PCs may then be rotated by applying a linear 
transformation according to the block labeled 170 to give rise to a new set of 
15 vectors called the rotated PCs (RPCs). Each of these RPCs contains all of the 
original quantities, and each vector can be defined by the largest absolute value 
coefficients. A new matrix [a] may be formed which contains the p rotated PC 
loadings for the reduced set of PCs. 

[0048] The PC model [z]=[F]/(| r can be rewritten to solve for the principal 
20 component scores [f] according to the block labeled 180, as follows: 

[0049] The p PC scores computed from the scaled data matrix and the 
25 RPCs became the predictors for the logistic regression, according to the block 
lebeled 190. The principal component scores [f] used in the logistic regression 
may be computed from the scaled data matrix [z] and the principal component 
loadings from the development matrix [a]. These scores may then be used in 
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the logistic regression to develop and test the weather event probability 
forecasts. 

[0050] In order to derive a probabilistic forecast for a binary event, (i.e., 
lightning either does or does not occur, or four inches of rain either does or 
does not occur), the invention may use logistic regression to fit the RPCs to 
a non-linear equation having an exponential in the denominator. Referring 
now to FIG. 6, this regression method produces an S-shaped curve bounded 
by 0 and 1 , which is suitable for results that are probabilities, as shown in 
graph 600. In addition to always producing a prediction between 0 and 1, 
the S-shaped curve of the logistic function, f(W), combines the contributions 
of the predictors and produces results that are slowly changing 610 for low 
values of w until a threshold is reached, then rises rapidly 620 for 
intermediate values of w and finally levels off 630 for high values of w. 
Historically, the parameters had to be estimated from an iterative technique 
that was computationally intensive, which limited its use in statistical weather 
forecasting. However, many statistical computer software packages now 
include logistic regression. 

[0051] Logistic regression has as input a binary predictand having the value 
of 1 whenever the predictand is true and 0 when it is false. The PC scores {S, 
to S p ) provided by selecting the PCs that contain most of the variance are the 

predictands. A generalized linear model may be used to derive the predictive 
equation having the following form: 

M 

where: 
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[0052] The maximum likelihood estimates of the coefficients (fi N ) are 

produced by an iteratively reweighted least squares approach. This common 
5 method is well known by those well versed in the field of statistics. 

[0053] The inventive method may be best illustrated and explained by 
example. An embodiment of the invention will be described for applying the 
method to a particular weather related event, namely, forecasting the probability 
of lightning strikes occurring. The inventive method will be used to develop a 

10 prediction system that produces three-hour grid point probability forecasts of 
cloud-to-ground lightning for the western United States. It should be 
understood that this illustrative embodiment does not serve as a limitation of 
the method but to illustrate and explain the general concepts that have been 
described previously. 

15 [0054] To fully appreciate the concepts and principles embodied by the 
invention as it is applied in this embodiment, a model may be established which 
describes the environment in which the method is applied. The features of this 
model may be dictated by the available data sets, the measurements that they 
represent, and the geographical areas to which they apply. It shall be 

20 necessary to define certain weather-related terms to better understand the 
particular data model and its assumptions. Some of these definitions and 
assumptions may be elemental to the field of meteorology, but are included for 
completeness. 

[0055] A sounding may be defined as a plot of the vertical profile of 
25 temperature and dew point (and often winds) above a fixed location. Soundings 
are used extensively in severe weather forecasting, e.g., to determine instability, 
locate temperature inversions, measure the strength of the capping inversion, 
obtain the convective temperature, etc. 



Page 17 of 48 



PATENT 
03-00076-03 



[0056] A cap (or capping inversion) may be defined as a layer of relatively 
warm air aloft (usually several thousand feet above the ground) which 
suppresses or delays the development of thunderstorms. Air parcels rising into 
this layer become cooler than the surrounding air, which inhibits their ability to 
5 rise further. The cap thus often prevents or delays thunderstorm development 
even in the presence of extreme atmospheric instability. However if the cap is 
removed or weakened, then explosive thunderstorm development can occur. A 
cap is an important ingredient in most severe thunderstorm episodes, as it 
serves to separate warm, moist air below and cooler, drier air above. With a 

10 cap, either the air below the cap can continue to warm and/or moisten or air 
above it can cool; either process increases potential instability. Without a cap, 
either process results in a rapid dissipation of the available instability, frequently 
before the severity of instability becomes large enough to support severe 
weather development. 

15 [0057] The Level of Free Convection (LFC) may be defined as the level at 
which an air parcel lifted from the surface becomes positively buoyant. The 
Equilibrium Level (EL) is defined, for a sounding, as the level above the LFC at 
which the temperature of a rising air parcel again equals the temperature of the 
environment, i.e. the pressure at which a rising saturated air parcel encounters 

20 negative buoyancy; above that pressure, the air parcel is cooler than its 
surroundings. The height of the EL is the height at which thunderstorm updrafts 
no longer accelerate upward. Thus, to a close approximation, it represents the 
height of expected (or ongoing) thunderstorm tops. However, strong updrafts 
will continue to rise past the EL before stopping, resulting in storm tops that are 

25 higher than the EL. This process sometimes can be seen visually as an 
overshooting top or anvil dome. The default units for EL are millibars, 
abbreviated as mb. 

[0058] Convective Available Potential Energy (CAPE) may be defined as a 
measure of the instability of the atmosphere, i.e. the amount of energy available 
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for convection. It may be further regarded as the amount of energy available to 
a surface air parcel that has reached its LFC. It is an integrated measure of an 
air parcel's buoyant energy in units of joules per kilogram (J/kg). The CAPE is 
also defined mathematically as an integrated value of the positive energy, 
5 observed on a sounding, which is available to an air parcel as it freely rises 
between the LFC and the EL. The CAPE is therefore non-zero only if an LFC 
exists. CAPE is directly related to the maximum potential vertical speed within 
an updraft; thus, higher values indicate greater potential for severe weather. 
Observed values in thunderstorm environments often may exceed 1 ,000 J/kg, 

10 and in extreme cases may exceed 5,000 J/kg. However, as with other indices 
or indicators, there are no threshold values above which severe weather 
becomes imminent. CAPE is represented on a sounding by the area (often 
called "positive area") enclosed between the environmental temperature profile 
and the path of a rising air parcel, over the layer within which the latter is 

15 warmer than the former. Another terms related to CAPE are Most Unstable 
CAPE (MUCAPE), i.e. CAPE based on the most unstable air parcel in the 
lowest 300 mb, and normalized CAPE (NCAPE), i.e., CAPE divided by the 
depth of the CAPE with units of acceleration (ms" 2 ). 

[0059] Convective Inhibition (CIN) may be defined as a measure of the 
20 amount of energy needed in order to initiate convection. Values of CIN typically 
reflect the strength of the cap. They are obtained on a sounding by computing 
the area enclosed between the environmental temperature profile and the path 
of a rising air parcel, over the layer within which the latter is cooler than the 
former (this area sometimes is called negative area.) 
25 [0060] The data set chosen to illustrate this embodiment of the invention 
may consist of the hourly three-dimensional data sets generated by the RUC2 
analysis as a basis for generating thunderstorm probability forecasts. The 
RUC2 has the advantage of rapidly incorporating all hourly surface data and 
any other off-time (i.e. non-radiosonde time) observations. Previous attempts to 
30 generate thunderstorm probabilities have had to rely on developmental data 
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sets initialized only from the normal analysis times (00 to 12 UTC), and then the 
thunderstorm forecast became valid 6 to 12 hours after the analysis. Since the 
RUC2 is able to capture important changes (especially at the surface) on an 
hourly basis, this can translate into hourly updated model forecast parameters, 
5 which are important to the accurate prediction of thunderstorms. 

[0061] The RUC2 model became operational in April 1998, and it marked the 
first time that a three-dimensional analysis was available on an hourly basis. 
The RUC2 data also contained the parameters necessary to evaluate the 
conditions needed for convection (pressure, temperature, moisture, and wind). 

10 The environmental model used for this embodiment may incorporate the 
physical processes that (1) initiate, (2) organize, (3) sustain, and/or (4) 
terminate convection by computing a large number of relevant diagnostic fields 
from the RUC2 data. In addition, a wide spectrum of diagnostic quantities may 
be calculated and evaluated at any desired level(s) in the atmosphere that might 

15 be physically related to thunderstorm genesis and maintenance. In addition to 
single level diagnostic quantities, other vertically averaged diagnostic quantities 
may be computed throughout the depth of the atmosphere. 
[0062] Previous attempts at predicting lightning strikes have not taken into 
account a number of issues. First, past studies have used data that is tied just 

20 to the surface or an arbitrary level, such as the 850 millibar (mb) level or 700 mb 
level. This truncation for data collection at only certain arbitrary levels does not 
consider the effect of conditions at other significant levels, which may be 
particularly pertinent when thunderstorms or other hazardous weather are being 
dealt with. 

25 [0063] Second, previous attempts, in an effort to simplify the calculations 
necessary and the observations to be taken, have not taken advantage of 
CAPE calculations. CAPE gives an indication of the instability of the 
atmosphere, so that the higher the CAPE value, the more unstable the 
atmosphere is. The RUC2 model provides values for CAPE on a regular basis. 
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Also previous attempts have analyzed only the surface based CAPE, rather 
than the total Most Unstable CAPE (MUCAPE), i.e. CAPE based on the most 
unstable air parcel in the lowest 300 mb. 

[0064] Third, the capping inversion has largely been ignored by prior efforts 
5 at lightning prediction, where the capping inversion is defined as a layer of 
relatively warm air aloft (usually several thousand feet above the ground) that 
suppresses or delays the development of thunderstorms. The presence of a 
strong capping inversion has a marked suppressing influence on thunderstorm 
development. Three diagnostic measures are used to identify and represent the 

10 strength of the capping inversion, i.e. the Convective Inhibition (CIN, as 
measured in J kg' 1 or mV 2 ), the Normalized CIN (CIN divided by the depth of 
the CIN and measured in ms 2 ), and the strength of the capping inversion 
(measured in degrees C). There has been little progress in bridging the 
temporal gap between observations and analyses and the first time period when 

15 model data is valid (usually 3 to 6 hours after the start of any model forecast). 
The MOS system represents a statistical system that cannot change as the 
model itself improves. This invention may employ a Perfect Prog approach and 
represent a statistical system not tied to any model; it can thus continue to be 
used as the model is improved. A reasonable expectation would be that as the 

20 model improves, so, in theory would the forecasts produced by this invention. A 
characteristic of the Perfect Prog methodology is that, as the forecast time 
period gets closer to the model run time, the forecasts should continue to 
improve. 

[0065] Forecasters have become more skilled in being able to predict the 
25 probability that a thunderstorm will become severe given that the thunderstorm 
develops. However, the most difficult part of the forecast often is to predict the 
probability of thunderstorms developing in the first place. 
[0066] Use of the environmental model described for this embodiment may 
have several important advantages. First, in deriving probabilities for every 
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three-hour period, different processes that may dominate during different time 
periods may be identified and highlighted. By deriving separate equations for 
each three-hour period, each of the equations may be more likely to capture this 
difference. For example, such differences might be daytime versus nocturnal 
5 storms or early afternoon development versus early evening decay of storms. 
Second, while some predictors (i.e. variables or data parameters) may be active 
through all eight of the three-hour periods, it is not reasonable to expect all 
terms to be acting in the same way throughout the course of any 24-hour 
period. Further, this environmental model may account for the influences of 
10 various terrain features in the area being studied by incorporating a set of time- 
dependent lightning climatological predictors developed especially for the model 
for each area of interest. It also may incorporate many new quantities that take 
into account the forcing mechanisms for thunderstorms throughout the depth of 
the atmosphere. 

15 [0067] The present embodiment illustrating the invention uses lightning 
climatological predictors that are determined for the eight contiguous three-hour 
periods in a day and are updated every five days by using an average derived 
for each five day period, or pentads. These predictors include the average 
number of flashes, relative frequency (or probability) of one or more flashes, 

20 average percentage of flashes over all grid boxes that occurred in each box, 
and average percentage of the three-hour interval that lightning is observed (i.e. 
the number of 15 minute intervals with lightning divided by 12). No current 
thunderstorm predictive system has a lightning climatology comprised of three- 
hour climatologies per five-day (pentad) period. 

25 [0068] The geographical area of interest may be divided up into rectangular 
grid boxes within which 3-dimensional climatological observations may be 
taken. For example, predictive data obtained from the western United States 
may be measured according to grid boxes that are 40 km. on a side. The 
choice of dividing the area into grid boxes and the size of the grid boxes may be 

30 dictated by the data and the ease by which it may be analyzed. Other choices 
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may used for developing a set of data points, such as, for example, using 
discrete locations at which measurements are taken, without departing from the 
scope of the invention. 

[0069] The present embodiment of the invention uses a large number of 
5 predictors computed from a multitude of model levels (every 25 mb in the 
vertical). For example, the level of the most unstable air parcel below 500 mb is 
included. Although the AVN/MRF MOS equations have recently used the 
Surface Based Convective Available Potential Energy (SBCAPE), no other 
model based predictive system has used MUCAPE and Most Unstable air 

10 parcel Convective Inhibition (MUCIN) together as a fundamental component of 
its system for thunderstorm prediction. For example, the current thunderstorm 
probability forecasts available to NWS forecasters from the ETA model are 
based only on the Lifted Index from the ETA and do not take into account the 
capping inversion or any other instability parameters. When CAPE and CIN 

15 measurements are used together, then areas of positive CAPE can be 
eliminated simply because the capping inversion (as measured by the CIN) is 
too strong for convection to develop. 

[0070] Total MUCAPE between 0° and -40°C. in the present embodiment is 
sub-divided into layers that are 5°C. thick based on the environmental 

20 temperature (e.g. 0° to -5°C. f etc.) In addition, when there is positive CAPE, the 
lowest and highest layers are computed from the Level of Free Convection 
(LFC) to 0°C. and from -40°C. to the Equilibrium Level (EL), respectively. The 
CAPE in these layers is defined as Layer CAPE (LCAPE). These layers may be 
used to determine any relationship between them and the presence of lightning 

25 and the number of flashes. Also, the relationship of Normalized CAPE 
(NCAPE) (i.e., CAPE divided by the depth of the CAPE with units of 
acceleration (ms" 2 )) to lightning and storms with high numbers of flashes may be 
examined. The Best CAPE (BCAPE) and CIN may then computed from the air 
parcel with the "best", or larger, CAPE from either of (1) the mean mixing ratio 

30 and potential temperature in the lowest 50 mb or (2) an air parcel above the 
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lowest 50 mb and below 500 mb. The BCAPE and CIN can give a better 
representation of the CAPE and CIN when there are unrealistically large surface 
moisture or super-adiabatic layers near the surface and can also capture the 
contribution of warm, moist air above a frontal surface or nocturnal inversion. 
5 [0071] Based upon the definitions given above and the availability of 
measurable and calculated data derived from the RUC2 data, Appendix A 
details a set of more than 200 predictors (parameters) of lightning-related 
weather phenomena that were derived from the RUC2 Model and used in the 
example illustrated by this illustrative embodiment. Other data and/or 
10 environmental models might have been used without departing from the scope 
of the invention or of this embodiment. Referring to FIG. 3, the data variables 
contained in the RUC2 model were classified for convenience of discussion 
according to the different vertical levels associated with each variable. Table 1 
gives the definitions of the acronyms used in FIG. 3. 

15 



Acronym 


Definition 


UL 


Upper Level support 


SR 


Wind Shear terms and/or storm 
Relative winds 


AK 


Averages of vertical motion 
(Adiabatic/Kinematic) 


TH 


Thermodynamics 


LC 


Lightning Climatology 


TH/LC 


best CAPE (Thermodynamics * 
Lightning Climatology 


SD 


Surface Data 


LL 


Lower-Level forcing 


MU 


Moist absolutely Unstable layer 


LP 


Lapse Rates 



Table 1. Explanation of Data Categories in FIG. 3. 



[0072] Referring to FIG. 4, a general flowchart is presented which illustrates 
20 the overall sequence of steps to be performed according to this embodiment. 
According to block 410, one input to the model consists of a set of gridded 
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historical weather data that can be obtained from any computer model having a 
known resolution according to time and space. For example, the RUC2 model 
has a resolution of 40 km for its grid. Another input to the model, according to 
block 420, is a set of observations of the data parameters according to the 
5 same resolution as the gridded weather data. The inputs provided in blocks 
410, 420 are used according to block 430 to derive a set of predictors believed 
to contribute to the selected weather event to be forecasted in a highly non- 
linear fashion. The variable to be predicted (as a probability) is converted 
according to block 440 into a binary predictand for forecasting purposes. For 

10 the present embodiment that forecasts the probability of lightning, one 
predictand may be whether or not one or more flashes occur, one predictand 
may be whether or not ten or more flashes occur, one predictand may be 
whether or not 100 or more flashes occur, and so forth. Given the desired 
predictands, a principal component analysis (PCA) is applied according to block 

15 450 to the set of predictors derived in block 430 to derive a new set of 
predictors having a much smaller number of members. This new set of 
predictors may be used to account for a majority of the variance present in the 
original set of predictors. A linear transformation (rotation) is also performed to 
transform the new predictors into a unique set of physically as well as 

20 meteorologically significant predictors related to actual atmospheric processes 
and/or variables. For example, one output of the rotation is a predictor made up 
of all original instability measures. Another is a measure of average low-level 
forcing and integrated vertical velocities. Finally according to block 460, a 
logistic regression is applied to the new (rotated) set of predictors derived from 

25 the PCA in block 450 to derive the probability of the selected weather event 
occurring; this probability may now be calculated given a new gridded 
(independent) data set. 

[0073] It is instructive to review some of the considerations that were made 
in determining the data model from which the predictive equations of the 
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present embodiment were derived. The data model for this embodiment was 
developed as follows. The RUC2 and lightning data were derived from 
observations made in June, July, and August of 1999. A lightning climatology 
was developed separately, using the years 1995 to 1998. See Attachment 1 for 
5 a complete list of the parameter data set. Quantities for the period 00 to 03 
UTC were used. 

[0074] Although lightning in the West can be very active into early 
September, the lightning climatologies confirm that the months of June, July and 
August normally represent the most active months for thunderstorms. Finally, in 
10 interpreting the forecast probabilities, it is realized that lightning over a three- 
hour period in a grid box that is 40 by 40 km is likely to represent multiple 
thunderstorms. The forecast probabilities do not attempt to represent individual 
thunderstorms; rather, they may represent a group of storms within each grid 
box. 

15 [0075] During the three months of June, July, and August1999, a total of 70 
days had both RUC2 and lightning data. Out of the 70 days, 54 days were used 
to develop the regression equations. Sixteen days were picked at random from 
contiguous pentads and withheld as independent data in order to test the 
equations. Since data from only one summer were available, the equations 

20 developed should be redeveloped when more data become available to improve 
the stability of the equations. 

[0076] The original data matrix [x] was composed of 208 candidate 
predictors, or parameters, in the columns, and grid point data for all 54 days in 
the rows (i.e., 2530 grid points for 54 days, or 136,620 observations). Initially, 
25 the first statistical analyses were attempted using all 136,620 grid points, even 
those without CAPE. However, since the purpose of the example was to search 
for physical relationships between groups of predictors and lightning, only the 
grid points where the MUCAPE was present (greater than 0) and extended from 
the LFC to at least the -15 to -20°C range were included in the final 
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developmental matrix [x], which had 89,115 rows by 205 columns. Many of 
the predictors (Attachment 1) are only defined in the presence of CAPE, 
otherwise, default values had to be applied (e.g., LFC, EL, etc.) The CAPE 
computed from the Most Unstable Parcel (MUP) was used to represent an 
5 "upper-bound" for the CAPE calculations. In reality, the actual CAPE values of 
a rising parcel could be far less than the MUCAPE, since the air at the Most 
Unstable Parcel Level (MUPL) may not be representative of the air going into 
the updraft. 

[0077] The layer CAPE from -15 to -20°C was found to be a good 

10 discriminator between cases with CAPE and lightning and cases with CAPE and 
no lightning. In other words, when any amount of CAPE was present, in over 98 
percent of the cases when lightning occurred, the CAPE extended to at least the 
layer from -15 to -20°C. Therefore, in order to concentrate on the cases where 
lightning was more likely (layer CAPE greater than zero from -15 to -20°C), the 

15 rows of data with no CAPE in this layer were removed from the matrix, leaving 
89,115 rows of data. This resulted in only 150 out of 13,523 grid boxes (1 
percent) with one or more lightning flashes being discarded. However, 47,505 
grid boxes (39 percent) without lightning were excluded from the analysis as a 
result of this condition. 

20 [0078] All 208 candidate predictors were used in the analysis and predictive 
equation development. However, out of the original 208 candidate predictors, 
205 candidate predictors were used in the PCA because the I, J, and K grid 
point locations were excluded, since these are fixed values with zero variance. 
The fundamental PCA method may be called an "R mode analysis" since an R 

25 mode is defined where the columns of the matrix are the predictors and station 
data (grid points in this case) are the rows. The PCA is termed a "Hybrid R 
mode" because a normal R mode analysis is for one time period and in this 
case 54 days are included in the rows. 
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[0079] A correlation matrix [r] (205 x 205) was computed from the scaled 
data matrix [z] (89,115 x 205). The eigenvalues [a] and eigenvectors [e] of 
the correlation matrix were then calculated. The individual elements, also called 
coefficients or PC loadings, for each of these 205 eigenvectors were computed. 
5 Mathematically, there are as many eigenvectors as there are original 
parameters, if there are no linearly dependent variables. A plot of the 
eigenvalues illustrates graphically how the number of eigenvectors to be 
retained can be determined. The graph 200 given in FIG. 2, sometimes called a 
scree graph, orders the eigenvalues in descending magnitude along the Y-axis 

10 221 corresponding to the PC index, along with its corresponding parameter, 
along the X-axis 211. Several methods exist to determine how many PC's 
should be retained, with the scree graph being one of these methods. A 
purpose of the scree graph is to separate the steeply sloping portion of the plot 
(containing the signal and some noise) from the portion where the slope is much 

15 less (containing primarily noise). As shown in FIG. 2, a rapid change in the 
slope appears in this case after the first ten eigenvalues 220. This would 
suggest that the ten parameters corresponding to the first ten PC indices 210 
along the X-axis 211 account for most of the variability in the data. The goal of 
data reduction would therefore be achieved by reducing the number of variables 

20 from 205 to 10. However, if the procedure were to stop at this point, there is the 
possibility that different physical processes could be included in the first PC 
(since it explains the maximum variance), making a physical interpretation 
misleading. Such unrotated solutions may exhibit four characteristics that 
hamper their utility to isolate individual modes of variation. These four 

25 characteristics are as follows: (1) the inaccurate portrayal of the physical 
relationships embedded within the data input matrix, (2) the PC's dependence 
on the shape of the domain, and/or (3) the size of a smaller subdomain, and (4) 
inability to physically interpret the results even with large sample sizes. 
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[0080] Since one of the goals of the method is the physical interpretation of 
the underlying order, an orthogonal (varimax) rotation, or linear transformation, 
may be applied to the first 10 loadings. Tests were also conducted applying the 
rotation to the first 12, 14, and 16 PC's. Minimal differences were noted in the 
5 results, as the first 10 rotated PC loadings were nearly identical in each of the 
cases. Increasing the number of terms that were rotated from 10 to 16 
accounted for minimal data compression (i.e. 58% of the original variance 
represented by the 10 terms versus 69% for 16 terms). Therefore, only ten 
terms were rotated since this was where a distinct "break" in the scree plot 
10 occurred and little additional variance was explained in terms 11 through 16. In 
addition to the varimax orthogonal rotation, two other transformations were 
attempted (oblique rotation (oblimin) and Procrustes) with no significant 
differences in the loadings. 

[0081] The rotation of the leftmost 10 PC's shown in graph 200 of FIG. 2 
1 5 may provide a new set of rotated vectors (RPCs). Each of these RPCs contains 
all of the original quantities and each vector can be defined by the largest 
(absolute value) coefficients. For the analysis of this embodiment, the 
coefficients with a value below 0.40 were not used because they contribute little 
information. The new RPCs and the percent of original variance explained are 
20 shown in Table 2. The scores computed from these RPC loadings were then 
used in the logistic regression. 



25 
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Nr. 


RPC 
Terms 


Amt of Var 
contained from 
original data 


Parameter 


1 


T 


13% 


CAPE 


2 


T 


10% 


Lapse rates/relative humidity 


3 


D 


6% 


Upper-level potential vorticity, saturated 
and moist geostrophic potential vorticity 
(above 400 mb), tropopause pressure, 
temperature, 500 mb height, height of 0, 
-10, and -20°C, WBZ 


4 


D 


5% 


Low level forcing and vertical velocity 
through LFC 


5 


D 


5% 


U wind components, shear and storm 
relative terms 


6 


D 


5% 


V wind components, shear and storm 
relative terms 


7 


D 


3% 


Mid-level saturated geostrophic potential 
vorticity (below 400 mb), mid and low level 
cyclonic circulation 


8 


T 


4% 


Convective Inhibition/Capping Inversion 


9 


T 


5% 


Best CAPE*LTG Climo, LTG Climo, 
potential instability 


10 


D 


3% 


Forcing from LCL to LFC, mid-level 
vertical velocity, upper-level divergence, 
pressure at most unstable parcel level 






59% total 





Table 2. Individual RPCs from Orthogonal (Linear) Rotation 



[0082] Referring to FIG. 5, a set of predictors is shown according to the 
5 taxonomy given in FIG. 3 that has been calculated according to the method 
given in FIG. 4 for the occurrence of lightning strikes. In FIG. 5, the ten 
predictors selected according to the inventive method have been divided into 
the ten groups and the percentage of variance explained by each predictor. 
The total variance explained by all ten predictors accounts for 59% of the total. 
10 [0083] The logistic regression for predicting the occurrence of lightning 
strikes according to this illustrative embodiment may have as input the binary 
predictand FONE, which is 1 for one or more flashes of lightning and 0 
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otherwise. The PC Scores computed from the PCA are the predictands and 
may then be calculated from the results of the previous process, where the 
developmental PC scores (£,.) correspond to the 10 RPCs shown in Table 2. 

The maximum likelihood estimates of the coefficients may be calculated by 

5 the iteratively reweighted least squares approach, and substituting in these 
coefficients derived from the logistic regression in the previous equation for W 
gives: 

W = -2.24 + 0.375, + 0.23S 2 + 0.28S 3 + 0.46S 4 - 0.28S 5 + 0.33S 6 
10 -0.015 7 +0.03S 8 +0.86S 9 +0.51S, 0 

[0084] As has been demonstrated the present invention provides 
advantageous techniques for a weather prediction system that will produce 
probability forecasts of at least an hour for selected weather events of interest. 

15 As will be appreciated by one of skill in the art, alternate embodiments of the 
present invention may be provided as methods, systems, or computer program 
products. Accordingly, the present invention may take the form of an entirely 
hardware embodiment, an entirely software embodiment or an embodiment 
combining software and hardware aspects. Furthermore, the present invention 

20 may take the form of a computer program product which is embodied on one or 
more computer-usable storage media (including, but not limited to, disk storage, 
CD-ROM, optical storage, and so forth) having computer-usable program code 
embodied therein. 

[0085] The present invention has been described with reference to flowchart 
25 illustrations and/or block diagrams of methods, apparatus (systems) and 
computer program products according to embodiments of the invention. It will 
be understood that each block of the flowchart illustrations and/or block 
diagrams, and combinations of blocks in the flowchart illustrations and/or block 
diagrams, can be implemented by computer program instructions. These 
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computer program instructions may be provided to a processor of a general 
purpose computer, special purpose computer, embedded processor or other 
programmable data processing apparatus to produce a machine, such that the 
instructions, which execute via the processor of the computer or other 
5 programmable data processing apparatus, create means for implementing the 
functions specified in the flowchart and/or block diagram block or blocks. 
[0086] These computer program instructions may also be stored in a 
computer-readable memory that can direct a computer or other programmable 
data processing apparatus to function in a particular manner, such that the 
10 instructions stored in the computer-readable memory produce an article of 
manufacture including instruction means which implement the function specified 
in the flowchart and/or block diagram block or blocks. 

[0087] The computer program instructions may also be loaded onto a 
computer or other programmable data processing apparatus to cause a series 

1 5 of operational steps to be performed on the computer or other programmable 
apparatus to produce a computer implemented process such that the 
instructions which execute on the computer or other programmable apparatus 
provide steps for implementing the functions specified in the flowchart and/or 
block diagram block or blocks. 

20 [0088] While the preferred embodiments of the present invention have been 
described, additional variations and modifications in those embodiments may 
occur to those skilled in the art once they learn of the basic inventive concepts. 
Therefore, it is intended that the appended claims shall be construed to include 
both the preferred embodiment and all such variations and modifications as fall 

25 within the spirit and scope of the invention. 
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APPENDIX A. Set of Lightning Predictors Derived from the RUC2 Model 



Var 
# 


GEMPAK 
Name 


Description 


Category 


1 


K 


K-grid point 


SD 


2 


I 


West to east grid point 


SD 


3 


J 


South to north grid point 


SD 


68 


LMSL 


Laplacian of mean sea level pressure 


SD 


69 


PTND 


Sfc pressure (altimeter) change (3 hr.) 


SD 


70 


LTND 


Laplacian of sfc pressure change (3 hr.) 


SD 


71 


LTHA 


Laplacian of surface Theta 


SD 


72 


PRES 


Pressure at surface (altimeter based) 


SD 


73 


ZMEG 


Terrain upslope/downslope 


SD 


74 


FRNT 


Surface frontogenesis 


SD 


75 


DQVC 


Surface Q vector convergence 


SD 


76 


PMSL 


Mean sea level pressure 


SD 


77 


RELH 


Surface humidity 


SD 


78 


MFCN 


Surface moisture flux convergence 


SD 


79 


WDIV 


Surface wind divergence 


SD 


80 


TWAD 


Surface Theta-W advection 


SD 


81 


RVRG 


Surface relative vorticity 


SD 


186 


LHGT 


Laplacian of surface height field 


SD 


187 


SMXR 


Surface mixing ratio 


SD 


198 


UWND 


Surface U wind component 


SD 


199 


VWND 


Surface V wind component 


SD 


4 


IOTF 


LTG average number flashes 


LC 


5 


ITFR 


LTG spatial relative frequency 


LC 


6 


IONE 


Itg relative frequency 


LC 


201 


ITSM 


LTG climo-% of 3hr pd Itg occurred 


LC 
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Var 
# 

TT 


GEMPAK 
Name 

I1QI 1 1 w 


D scription 


L»aiegory 


204 


BOTF 


BECP*IOTF 


TH/LC 


205 

C-\J\J 


BTFR 

1— S II 1 \ 


BECP*ITFR 

VJ 1— >■/ 1 III IX 


TH/LC 


206 


BONE 


BECP*I0NE 


TH/LC 


207 


BTSM 


BECP*ITSM 


TH/LC 


46 


FCRH 
■ x> i \i i 


AVG humiditv from LCL to LFC 

/ \ V ^mj 1 IUI 1 II VI 1 V ¥ II W III Imm Ib * W b 1 


LL 


64 


RHAV 

i xi in v 


AVG humiditv from MUPL to LCL 

# % v \mf iiuiiiiwJty iiviii iviwi i— iw uvu 


LL 


105 


ASRH 

nui \i i 


AVG sub-cloud humiditv 

# \ V \m* WUU VlWUVl 1 IUI 1 1 IU 1 I J 


LL 


10Q 


MCAV 


AVG moist flux conva (MU PL-LCD 

/% V 1 1 IUIOI IIU/\ V/WlIVy yiVIWI mm, mm.\m* mm, J 


LL 


110 


WCAV 

V V \m*i\ V 


AVG wind diveraence (MUPL-LCL) 

#» V \m+ Wlllw. UIVvlUvl IWW ^IVIWI 1— m—\m*mm. J 


LL 


111 

III 


TCAV 


AVG Theta-W advection (MUPL-LCL) 


LL 


112 

1 I £— 


RCAV 


AVG aeostroohic rel vort (MUPL-LCL) 

/ \ V X^ U w WW 11 X/L/ 1 1 1 w 1 wl V VI 1 1 V 1 \m? 1 k_ L— V—/ L— J 


LL 


113 


FCAV 


AVG frontoaenesis (MUPL-LCL) 

f \ V \^ 1 1 VI 1 IwMVl 1 WWl W 11 VI 1 mm, mm, mm, t 


LL 


114 


LCAV 


AVG Theta LaDlacian (MUPL-LCL) 


LL 


115 


HCAV 


AVG heiaht (pressure sfc) (MUPL-LCL) 


LL 


116 


QCAV 


AVG Q vect converaence (MUPL-LCL) 

§ \ w VK ▼ V* VI 1 1 w W I W 1 1 W 1 I w 1 ■ mm, mm, mm, § 


LL 


117 

i i r 


MFAV 

1 VI 1 /» V 


AVG moist flux conva (LCL-LFC) 

/ \ W >mf III I W * II W. 1 l» 1 mm, m^ mmm I g 


LL 


118 


WFAV 

¥ V 1 #l V 


AVG wind diveraence (LCL-LFC) 

1 \ V V V 1 1 1 Vrfl VI 1 V W 1 W 1 1 WW 1 tmm ^mT mm, mm, I f 


LL 


119 


TFAV 

1 1 / V V 


AVG theta-W advection fLCL-LFC) 


LL 


120 


RFAV 


AVG aeostroohic rel vort (LCL-LFC) 

rVV V UuUvllV/li/lllv 1 vl VUl l ILvL l—l w J 


LL 


121 

1 Cm 1 


FFAV 

\ 1 l\ V 


A\/n frnntonenp^i^ (LCL-LFC) 

#\ v y II Ul IIU^CI ICOIO ^LVL l—l 


LL 


192 


1 FAX/ 


AVG Theta Lanlacian (LCL-LFC) 

nv w 1 1 Ida L.CipiOwlCll 1 yLUL l—l w f 


LL 


123 


HFAV 


AVG height Laplacian (LCL-LFC) 


LL 


124 


QFAV 


AVG Q vect convergence (LCL-LFC) 


LL 


132 


AWLC 


AVG ageostrophic wind (MUPL-LCL) 


LL 


133 


AWLF 


AVG ageostrophic wind (LCL to LFC) 


LL 


134 


VUAV 


Change-ageostrohic shear (LCL-400 mb) 


LL 


135 


GUAV 


Gradient vertical motion (LCL-400 mb) 


LL 
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136 


VLAV 


Chanae-aaeostroDhic shear (LCL-LFC^ 


LL 


137 


GLAV 


Gradient vertical motion (LCL to LFC) 

1 UU IW 1 1 W V V* 1 hi WM 1 III \S wl V/ III Va IV I 1 


LL 


188 


MXTE 


Theta-E at MUPL 

III V* V_ Wl k 1 V 1 V/ 1 


I L 


7 


MLIM 


Most neoative LI (MUP) from anv level 


TH 
■ ii 


8 


MLIO 


Pressure level of MLIM 

I I V^WVf VI 1 V ■ V V W 1 \S 1 1 V 1 k* I 1 V 1 


TH 
■ ii 


9 


MLIH 


Heiaht of MLIM 

1 1 V# 1 1 l\ V# 1 IVIL.IIVI 


TH 
iii 


11 


MELP 


Pressure of the EL usina the MUP 


TH 
iii 


12 


MELH 


Heiaht of the EL usina the MUP 

1 1 V/IVjl 11 VI VI IV WVII VJ * IV III VI 


TH 
■ ii 


13 


MELT 


Temperature (°C) of the EL 


TH 


14 


CDPM 


Cloud depth (meters) from LCL to EL 


TH 


15 


FCDP 


Depth LCL to LFC (when LFC present) 


TH 


16 


MLPP 


Originating pressure level of MUP 


TH 


17 


MBST 


Lift (mb) for MUP to reach saturation 


TH 


18 


MMPP 


Pressure (mb) Max air Parcel Level (MPL) 


TH 


19 


MMPH 


Heiaht (meters) of MPL 

1 1 1 1 II 11 1 1 1 f V<# 1 1 V 1 1 w 


TH 
iii 


20 


NCAP 


Normalized Most Unstable CAPE (m s" 2 ) 


TH 
■ ■I 


21 


BECP 


Best CAPE-best-50 mb ava or level above 

fc^ V* V L V«r # %l ^— V/vvV vw 1 1 IU V VI VI IvVvl V4vv/Vv 


TH 
■ ■I 


22 


BECN 


Convective Inhibition (CIN) for BECP 

vwi ivvviivv 1 1 ■ I 1 1 I ii v/i I i VII 11 1 vl I— » k— w i 


TH 
1 1 1 


23 


EMCP 


Elevated most unstable CAPE above cao 

VjI v V UlvVi 1 1 1 vvl Wl 1 IvlvlMlV Vni 4— V4 v v V v UQk/ 


TH 
iii 


24 


EMCN 


CIN for EMCP 

V/ll 1 1 vl VjIVIVI 


TH 
iii 


25 


UUCP l 


Unstable and uncaDoed CAPE 

V*l IWlUVIv VI 1 1 VI VII IVUUpvU W/ \l 1— 


TH 
iii 


26 


MUCP 


Most Unstable CAPE (MUCAPE) 


TH 


27 


MUCN 


CIN-most unstable CAPE (MUCAPE) 


TH 


28 


MULI 


LI for most unstable air parcel at 500 mb 


TH 


29 


MLFO 


Layer CAPE (LFC to 0°C) 


TH 


30 


MM05 


Layer CAPE (0 to -5°C) 


TH 


31 


MM10 


Layer CAPE (-5 to -10°C) 


TH 
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32 


MM15 


Laver CaDe (-10 to -15°C) 


TH 


33 

wv 


MM20 


Laver CAPE (-15 to -20°C) 


TH 


34 


MM25 


Laver CAPE (-20 to -25°C) 


TH 


35 


MM30 

1 V 1 1 V 1 ww 


Laver CAPE (-25 to -30°C) 


TH 


36 


MM35 

1 VI 1 VI ww 


Laver CAPE (-30 to -35°C) 

LuVwl 11 U. 1 ww I w ww 1 


TH 


37 


MM40 

i vi i vi~ v/ 


Laver CAPE (-35 to -40°C} 

LQy vl ^// \i Li I W W Iw ? w 1 


TH 


38 


MM4E 

IVI 1 Vl*t l_ 


Laver CAPE f-40°C to EL^ 

L_CI V wl \_//\l I— I ~ w lu 1— L- / 


TH 

ill 


39 


MLCT 

1 V 1 1_ \_y 1 


Cloud ha^p iMCL ^ tpmnpraturp (°C^ 


TH 
1 1 1 


40 


MH10 

IVI 1 1 1 w 


Heiaht fmeters-MSL) of-10°C 

1 IvIM III III IWlWl w IVI ^— I \J 1 1 \J w 


TH 


43 


MH20 

IVII 1 t— w 


Heiaht fmeters-MSL) of -20°C 

1 1 V./ 1 VJ Ilk till | W IVI W k_ 1 \J 1 W \^ 


TH 


50 


MUCF 


Most unstable CAPE from LCL to LFC 

1 V 1 WW k Vi 1 IWVUmIV \ 1 ft— 1 1 \J% 1 1 ftaa V./ 1— IV ft— 1 


TH 


53 


MCAP 


Caooino inversion for MUP (°C) 

vM 1^ 1^ III \A 1 1 1 v 1 l ll l x^ l IVI 1 \ ' J 


TH 


54 


LMCN 


Laolacian of CIN 


TH 


55 

V/ W 


EFDF 


Difff meters) EL to freezina level 


TH 


56 


EWDF 


DifWmeters) EL to wet-bulb zero 

kflf 1 1 1 1 1 1 1 V 1 \^ J mmmm %mm V W V V X<* W V* ■ fc* X^ I X^ 


TH 


57 


WBZH 


Heiaht of wet-bulb zero (meters-MSL) 

1 1 k^ 1 Vj Ilk 1 If k »«* W 1 fc* «— 1 11 ■ 1 k^ fc>^l W 1 » 1 f 


TH 


58 


FZLH 


Heiaht of freezina level (meters-MSL) 

1 1 ■ y Ilk 1 1 1 1 1 1 ■ \^ V III 1 wlwl W IVI Iw / 


TH 


59 

w w 


DNCP 


Downdraft CAPE (Jlko) 


TH 

ill 


60 


TEDF 


Thpfa-F riiffprpnpp mid-lnpx/pl^ tn ^fr 

1 1 1 w LC4 ft— 1 1 1 wJ vl 1 ww 1 1 IIVJ Iww V wlO Iv wl w 


TH 
1 1 1 


61 


MXMX 


Mivinn ratin for MUP 


TH 


62 


OTRN 


MXMX*storm rel inflow at MUPL 

IVI/\IVI/\ OiVJIIII 1 wl IIIIIW/W Ol IVlvyi L 


TH 


63 


XTRN 


MXMX*wind speed at MUPL 

1 f 1/ \ l ff 1/ X Iff II 1 \*A \*S V<J \A IVI 1 ka 


TH 


82 


MLCH 


Cloud base height (meters-msl) for MUP 


TH 


83 


MLFH 


Height of LFC for MUP 


TH 


106 


INPW 


Inches of precipitable water 


TH 


189 


MULP 


Laplacian of moist unstable CAPE 


TH 


190 


THEL 


Laplacian of Theta-E from MUPL 


TH 
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191 


TRPK 


TemD (°K\ at troDODause ("2 PV" units ^ 

1 CI 1 1 U 1 IXJ CI l 11 Wk/WUUUww I £- 1 V Ul 1 1 to J 


TH 
iii 


192 


TRPT 


l^pntrnnir ^nrfarp of trnnnn^iiQp 

iociiliuuiu oui ictuc ^ f\j ui ii upupauoc 


TH 

i r i 


195 

1 C/ U 


MCNP 


Prp^urp IpvpI of ranninn inx/pr^inn 


TH ' 
1 1 1 


196 


ML CP 

IVI Lwl 


Prpcciirp IpvpI of 1 CI for MUP 

1 1 wwOUl v ICVUl UI 1— \S 1— IUI 1 V 1 U 1 


TH 
1 1 1 


197 
i u i 


ML FP 

ivi Li i 


Prp^urp IpvpI nf LFC for MUP 

i icoouiu icvci ui i— 1 w i ui ivi ur 


TH 
mi 


200 


NCIN 


Normalizpd CIN f'm/s 2 ^ 


TH 
1 1 1 


202 


CTOT 


Laver CAPE from LFC to -20°C 


TH 
iii 


203 


CZ20 


Laver CAPE from 0 to — 20°C 

uayui w/»i i— ii will u wj v*/ 


TH 
■ 1 1 


208 


BMCR 


Ratio BECP to MUCP (0 if MUCP=0^ 

i \aiiu WLWI LU IVIUVI IV II IVI wVI Uf 


TH 
1 1 1 


41 

"T 1 


ELSU 


U shpar romnonpnt from FL to LCL 

W ul IvQI vUl 1 IL/V/I lv/1 11 IIUI II I— L LU LvL 


SR 

vvl X 


42 


ELSV 


V shear comoonent from EL to LCL 

V wl IwUl IA/I 1 IUUI Ivl 11 II Will 1— I— IV/ L_ V_/ L_ 


SR 1 

vl X 


44 


FCSU 


U shear comoonent from LCL to LFC 

\-J wl IvUI VV/I 1 1 VJ \J 1 Ivl III II vl 1 1 ^V./w IW 1 v»/ 


SR 

VVI X 


45 


FCSV 

1 WW V 


V shear comoonent from LCL to LFC 

V 0 1 IV/CII wl 1 IUUI Ivl 1 1 IIUI II LvL LU LI W 


SR 

WIX 


84 
u~ 


UPMW 

W 1 IVI V V 


Sfo to fi km nrp^^urp-wpinhtprl I J romn 

v lu u rx. iii [jicoouic vvciyi lieu \j uui 1 iu 


SR 

wrx 


85 


VPMW 

V 1 IVI V V 


Sfr to fi km nrp^<5iirp-wpiohtpri V romn 

vlv LU U f\l 1 1 Ul COvUI v VVUIUI llVv V vUl 1 IU 


SR 

urv 


86 


UCBL 


Clond-hparinn lavpr LJ mpan wind romn 

viuuu ucgii ii ly lay ci w iiicciii vviiiu uuiiiu 


SR 

Ol X 


87 


VCBL 

V VvLV L 


ClouH-hparinn lavpr V mpan wind romn 

vviuuu ucai ii ly lay vi v 1 1 icai 1 wniu uuiiiu 


SR 

vvl X 


88 


U3SV 


Surfarp to 3 km 1 J Qhpar romnonpnt 

UUI luvU LU U i\l 1 1 w Ol ICO 1 UUI 1 IUUI ICI IL 


SR 

VVIX 


89 


V3SV 

V Uw V 


fiurfarp to 3 km V ^hpar romnonpnt 

UUI luw LU U l\l 1 1 V Ol IvQI UUI 1 IUUI ICI 1 L 


SR 

VVlX 


90 


U6SV 

WWW V 


fiurfarp to fi km I J ^hpar romnonpnt 

UUI law LU U l\l 1 1 W Ol ICCII UUI 1 IUUI ICI 1 L 


SR 

WIX 


91 

V 1 


V6SV 

V U W V 


Surfarp to fi km V ^hpar romnonpnt 

vUl luvC LU U r\l 1 1 V Ol ICCII UUI 1 1 UUI ICI IL 


SR 

wrx 


92 


ULSV 


Low level to 6 km U shear comoonent 

^— v W w IV W V 1 VV \S 1X1 1 I Wl 1 Vvl 1 VX^I 1 1 VV ■ IV lit 


SR 


93 


VLSV 


Low level to 6 km V shear component 


SR 


94 


BRNS 


Bulk Richardson Number Shear (m 2 /s 2 ) 


SR 


95 


USSR 


Surface to 2 km storm rel. U wind comp 


SR 


96 


VSSR 


Surface to 2 km storm rel. V wind comp 


SR 


97 


U4SR 


4 to 6 km storm rel U wind comp 


SR 
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98 


V4SR 


4 to 6 km storm rel V wind comp 


SR 


99 


U9SR 


9 to 1 1 km storm rel U wind comp 


SR 


100 


V9SR 


9 to 1 1 km storm rel V wind comp 


SR 


101 


UASR 


Anvil level storm rel U wind comp 


SR 


102 


VASR 


Anvil level storm rel V wind comp 


SR 


103 


UMXP 


U wind comDonent at MUPL 


SR 


104 


VMXP 


V wind component at MUPL 


SR 


125 


INFL 


Storm relative inflow (maanitude) 

■ V ■ 9 9 M* P%M mm w Mr ■ ■ ■ ■ ■ WW 111 1 MM ■ I ■ Vw \Jl \4 m 


SR 


126 


OTFL 


Storm relative outflow (magnitude) 


SR 


179 


SH38 


300 to 850 mb speed shear 


SR 


48 


VKCL 


Kinematic vertical velocity at LFC 


AK 


49 


WCL 


Adiabatic vertical velocity at LCL 


AK 


51 


VKFL 


Kinematic vertical velocity at LFC 


AK 


52 


WFL 


Adiabatic vertical velocity at LFC 


AK 


108 


WAV 


AVG adiabatic vert vel (LCL-400 mb) 


AK 


127 


WLC 


AVG adiabatic vert vel (MUPL-LCL) 


AK 


128 


VKLC 


AVG kinematic vert vel (MUPL-LCL) 


AK 


129 


WLF 


AFG adiabatic vert vel (LCL-LFC) 


AK 


130 


VKLF 


AVG kinematic vert vel (LCL-LFC) 


AK 


131 


VKAV 


AVG kinematic vert vel (LCL-400 mb^ 


AK 

/VIA 


66 


TELP 


Moist Absolutelv Unstable Laver 


MU 

IVI w 


67 


MMBD 


Depth of MAUL (mb) 


MU 


47 


FCTL 


Theta difference MUPL to LFC 


LR 


65 


THLA 


Theta difference from MUPL to LCL 


LR 


107 


LLLR 


Lower level lapse rate-sfc to 3 km 


LR 


138 


LP78 


Theta-change/press diff 700 to 850 mb 


LR 


139 


LC78 


LP78*AVG converg 700 to 850 mb 


LR 
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140 


Ql 78 
OLf O 


Qat ThAta.F/nrpQQ Hiff 700 to 850 mh 

OCl l II lUld U/ Ul UOd III 1 UU LU UUU WXVJ 


LR 


141 
It I 


Ml 7ft 


Th^tfl-F/nrPQQiirp rliff 700 to 850 mh 
1 1 mid c/ui cooui c uni # iu uju iiiu 


LR 


144 


1 Pft7 


Thpfa-r^hannp/nrPCQ Hiff 800 tn 700 mh 
1 MCLd ~ui idi iyc/ ui coo uiii uuu iu iuu iiiu 


LR 


1 AR 
1 40 


1 Pft7 
LOOf 


1 PR7*A\/^ rnnvprn R00 tn 700 mh 
Li ur nVO ouiivciy uuu iu i uu iiiu 


LR 


140 


Ql A7 


Qat Thofo.F/nroQQ Hiff 800-700 mh 
Odl 1 ilCld-c/|Jlcoo uiii uuu- iuu iiiu 


LR 

l—l X 


147 
14/ 


Ml 87 
IVI LD f 


ThptsuF/nrPQQ Hiff R00 tn 700 mh 

1 1 It? Id C/ Ul COO UIII UUU WJ 1 \J\J IIIU 


LR 


159 


1 P58 


Thpfa-rhannp/nrpQ^ Hiff 500 tn R00 mh 

1 1 ICld Ul idi ly c/ yji coo uin uuu iu iiiu 


LR 


15^ 


1 P58 


1 P58*AVfi rorwpro 500 to 600 mb 

i_i ju nvu uuiivoiy jvu iu iiiu 


LR 


1^4 

I U*T 


Ql 5fi 


Qflt Thpta-F/nrp<5<; Hiff 500-fiOO mh 

Odl 1 I icia L./|jicoo uiii Jwu 1 1 iu 


LR 


1RR 

I uu 


Ml 58 


Thpffl-F/nrp<;<; Hiff 500 tn R00 mh 

1 lit? Id L-/UICOO UIII UUU IU UUU IIIU 


LR 

l— 1 X 


I uo 


Tl 75 


Tpmnpratiirp Hiff-700 to 500 mb (°C) 

1 CI 1 IUCI dlUI C UIII / UU IU JUU MIL/ ^ vy 


LR 


1 RQ 


1 PQ4 


Thpfa-nhannp/nrPQQ Hiff 400 mh to ^fc* 
1 1 icid oi idi lyc/ ui coo uiii *tuu iiiu iu oiu 


LR 

1— 1 X 


1fiO 


1 P48 


Thpta-rhannp/nrpQQ Hiff 400 to 850 mh 

1 1 icid ui idi iyc/ Ul COO UIII tuu IU UUU IIIU 


LR 

^>l X 


181 


1 P47 


Thptsurhannp/nrpQQ Hiff 400 to 700 mh 
1 1 icid ui idi lyc/ui coo unt ■tuu iu r uu iiiu 


LR 

LI X 




1 PAR 


Thpfa-nhannp/nrPQQ Hiff 400 to 500 mh 
1 1 icid ui idi iyc/ ui coo uiii *tuu iu uuu iiiu 


LR 

L_ 1 X 


18^ 
I DO 


1 P47 


1 P58*AVf5 ronvpro 500 to R00 mh 
L-i uu nvu uuiivciy uuu iu uuu iiil/ 


LR 

L.I X 


AAA 
I D4 


1 P45 


1 P45*AW^ ronx/prn 400 to 500 mh 

Lr*tJ nvu uuiivciy *tuu IU UUU IIIU 


LR 

LI X 


I DO 


Ql Q4 


Qat Thpta-F/nrpQQ Hiff 400 mh to Qfp 

Odl 1 l ICld L_7 Ul coo UIII *+uu IIIU IU OIU 


LR 


1£A 
I DD 


Ql Aft 


Qat Thpta-F/nrPQQ Hiff 400 to 850 mh 

Odl 1 IIC Id "L-/UICOO UIII *tUU iu uuu iiiu 


LR 

L.I X 


187 


Ql 47 


Q«t Thptfl-F/nrp<;^ Hiff 400 to 700 mh 
Odl 1 1 icid c/ui coo uiii tuu iu i uu iiiu 


LR 

L_l X 


I DO 


Ql 45 


Qat Thpta-F/nrPQ^ Hiff 400 to 500 mh 

Odl 1 1 ICld L-/ Ul coo uiii tuu iu uuu iiiu 


LR 

L«l X 


169 


MLS4 

IVI LiVr~ 


Theta-E/oress diff 400 mb to sfc 


LR 


170 


ML48 


Theta-E/press diff 400 mb to 800 mb 


LR 


171 


ML47 


Theta-E/press diff 400 mb to 700 mb 


LR 


172 


ML45 


Theta-E/press diff 400 mb to 500 mb 


LR 


173 


SP45 


Sat Geostrphic PV-400 to 500 mb 


LR 


174 


MP45 


Moist Geostrphic PV-400 to 500 mb 


LR 
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175 


SP34 


Sat Geostrophic PV-300 to 400 mb 


LR 


176 


MP34 


Moist Geostrophic PV-300 to 400 mb 


LR 


10 


TRPP 


TroDODause (as defined bv 2 PV units) 


UL 


142 


SP78 


Sat Geostroohic PV-700 to 850 mb 


UL 


143 


MP78 


Moist Geostroohic PV-700 to 850 mb 


UL 


148 


SP67 


Sat Geostroohic PV-600 to 700 mb 


UL 


149 


MP67 


Moist Geostroohic PV-600 to 700 mb 


UL 


150 


H5HT 


500 mb heiaht field f meters } 

w X* til fcy 1 1 %y 1 %J 1 11 II >^ 1 V* till \* vv 1 w J 


UL 


151 


H5LP 


Laolacian of 500 mb heiaht field 

^m\**S 1 UWIV4 ■ 1 I V*^ V/ V III 1 1 VS 1 lib II 1 \Ji 


UL 


156 


SP56 


Sat Geostrophic PV-500 to 600 mb 


UL 


157 


MP56 


Moist Geostrophic PV-500 to 600 mb 


UL 


177 


UREL 


300 mb U wind component 


UL 


178 


VREL 


300 mb V wind component 


UL 


180 


AUDV 


AVG upr-lvl divergence (250 to 400 mb) 


UL 


181 


AUAW 


Mag ageostrophic wind (250 to 400 mb) 


UL 


182 


PVOR 


Potential vorticity (PV) (250 to 400 mb) 


UL 


183 


PVAV 


PV advection by 300 mb wind 


UL 


184 


SP23 


Sat Geostrophic PV 200 to 300 mb 


UL 


185 


MP23 


Moist Geostrophic PV 200 to 300 mb 


UL 


193 


TRPU 


U wind component at tropopause level ! 


UL 


194 


TRPV 


V wind component at tropopause level 


UL 
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